Entanglement invariant for the double Jaynes-Cummings model 



3 



> 



INTRODUCTION 
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(Dated: October 10, 2008) 

We study entanglement dynamics between four qubits interacting through two isolated Jaynes- 
Cummings Hamiltonians, via an entanglement measure based on the wedge product. We compare 
the results with similar results obtained using bipartite concurrence resulting in what is referred 
to as "entanglement sudden death". We find a natural entanglement invariant under evolution 
demonstrating that entanglement spreads out over all of the system's degrees of freedom that be- 
come entangled through the interaction. We also provide an analysis why certain initial states 
loose all their entanglement in a finite time although their excitation and coherence only vanishes 
. asymptotically with time. 
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■ Entanglement plays a key role in quantum information processes [1| and therefore it is important to study entangle- 
ment dynamics in different scenarios. The simplest situation, two-qubit entanglement dynamics, has been extensively 
1 ■ studied in different contexts @, S, H 1, & 0, 1, ©, M, EL El ■ 

When interacting with a reservoir, one would naively expect the entanglement between the two qubits to vanish 
. asymptotically. However, for certain initial entangled states, the entanglement can vanish completely in a finite time. 
Q_i' This is often referred to as "entanglement sudden death," see for example [!, 0, H, @, 0, H, [1[ and the references therein. 
+L» . A recent experimental demonstration was presented in [l(| and an open systems analysis was done in [H, 0, H, IS 0] • 
Nevertheless, if we include the "reservoirs" in the studied system and consider the full entanglement between two non- 
interacting partitions of the system, one would expect the entanglement to be preserved, and therefore an associated 
entanglement invariant should exist. In addition, considering the reservoir as part of closed system until the moment it 
is ignored (and in mathematical terms "traced over" ) , one can reap insights into both the qualitative and quantitative 
transfer of excitation (and associated entanglement) from the atoms to the reservoir. 

An important bipartite interaction is described by the Jaynes-Cummings (JC) model [l3| that describes, in a concise 
and elegant way, the near-resonant interaction between a single two-level atom and a single-mode quantized field. If 
the initial atom-field system contains a single excitation, the system is a model of a two-qubit system. Since the JC 
qq , model is excitation number preserving, the system will always stay within the two qubit Hilbert space (but of course 
£f) • such a JC model spans only the one-excitation subspace of the full two-qubit space). However, since this model is one 
, -J of the few exactly solvable models in quantum physics it has been exploited for studying the dynamics of entanglement 
0- 

Recently, a double JC model has been proposed in this context [1,0]. The model consists of two separate JC- model 
systems (atom A interacting only with the cavity field a and similar for the atom B and the field b), where it has 
been assumed that the systems are identical. (Note that this model is applicable to any one-excitation, two-qubit 
• i-H , system that is linearly coupled, e.g., to the experiment in flpj where fields couple pairwise to each other rather than 
atoms to fields.) A major reason this particular interaction has been chosen is because it is local to subsystems Aa 
and Bb. If the atoms, or the fields, couple to each other, the coupling will alter the entanglement between the system 
partitions Aa and Bb in general, and subsequently, if the fields are traced over, between A and B. The whole point 
with these studies, however, is to study the entanglement dynamics between A and B in absence of any coupling, 
direct or indirect, between them. 

The focus of interest has been the pairwise entanglement dynamics in terms of the concurrence [l5[ between the 
initially entangled atoms. Through the JC interaction they may become unentangled through the excitation transfer 
to the initially unexcited fields which are traced over after the interaction. In Ref. [sj, in particular, the authors 
study entanglement between the two atoms and they find that for the initial state |^(0)) = cos a |||) + sin a |J,j), 
where we have used the notation \\) A ® \i) B = \M), etc. and |f) (||)) denotes the atom's excited (ground) state, 
the concurrence Cab behave in a harmonic oscillatory manner. Translating this into a dissipation language, the 
entanglement vanishes asymptotically with increasing coupling to the reservoir. In contrast, the concurrence Cab of 
the state |<£(0)) = cos a |f |)+sina | J, J.) for certain values of a, specifically | tanaj < 1, falls rapidly and non-sinusoidally 
to zero in a finite time and remains zero for some time. In a dissipation language this means that the entanglement 
will vanish in a finite time although both the atomic excitation and the atomic coherence decay asymptotically. 

In order to study the transfer of entanglement between the atom and the reservoir the authors extend the work in 
Q and study all the 6 concurrences, Cab, Cao,, C B b, C a b, CAb, and Cs a , for the four-qubit system in Q. They find that 
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for the state |-0(O)) the sum of the concurrences between the atoms and the fields, Cab + C a b, is constant under the 
JC evolution. This is not so for the state |</>(0)}, but another function of the six pairwise concurrences and the initial 
state (parameterized by a), Cab + C a b + {Cau + CBb)\ tana| — (Ca& + Cbo), is an entanglement invariant. 

In [l2| the dynamics for the initial state IV'(O)} in an equivalent model, two separately systems composed by two 
two-level systems in a dipolc-like interaction, is considered. Here, the authors consider the effect of different coupling 
constants in the separately systems. However, the pairwise entanglement dynamics is expressed in terms of the 
negativity [l6| and its relation to with the energy transfer. 

Motivated by these studies, and in the hope that a more general entanglement invariant can be found to illuminate 
the transfer of entanglement to the system's different parts, we study the same four-qubit system, where we do 
not assume that the systems are identical. To be able to consider all possible bipartite entanglement we use an 
entanglement measure introduced by Heydari in [13] which is based on the wedge product. We take all the different 
possible partitions of the four-qubit system into account, and we find an entanglement invariant which does not depend 
on the system parameters or on the initial state, provided that it belongs to a class of pure states denoted "X" -states 
in including both the state |-0(O)} and |^>(0)}. The invariant shows that the fields become entangled with all other 
parts of the system, all of which is "destroyed" (or rather ignored) when treating the fields as reservoirs. 



II. THE MODEL 



Consider a model consisting of two two-level atoms A, B, each interacting with a single-mode near-resonant cavity 
field denoted a and b, respectively. Following HQ, we will assume that each atom-cavity system is isolated and that 
the cavities are initially in the unexcited state while the atoms are initially in an entangled state. The dynamics of 
this model is given by the double JC Hamiltonian 

H tot =H A + H B , (1) 
where the Hamiltonians (under the rotating wave approximation and setting h = 1) are [l3j . 

H k = Vk{a\a k + 1/2) + ^-a\ + g k (ala k _ + a k a k + ), (2) 

where k = A, B (where the letter case is to be interpreted as appropriate), v k is the field frequency, iOk is the transition 
frequency between the atomic excited and ground states, and g k is the coupling constant between the cavity field 
and the atom. The field annihilation operators are fife, and a± are the spin-flip operators defined by a_ ||) fe = \i) k , 
°- \l)k = °+ lT)fe = °! °+ IDj. = |T) fe) and (7^ is the atomic inversion operator, viz., u\ \\) k = ||) fc and a k z \\) k = - |J.) fc . 
As mentioned in the Introduction, if we have at most one excitation in each atom-cavity system, each such system 
will stay within a two-qubit space. Hence, since the two atom-cavity systems don't interact, the double JC model will 
result in a four qubit state (but again, spanning only a subspace of the whole four-qubit Hilbert space). 
The corresponding evolution operator for the Hamiltonian ^ is 

Uk = <r it&l {cos(f! fe i) (3) 



-it sinc(Qfct) 



where H k — Vk[a\a k + (1 + a k )/2)] is a constant of motion, proportional to the total number of excitations of system 
k, Afc = u>k — Vk is the detuning between the atom and the cavity for each system, sinc(a;) = x^ 1 sin(x), and 

1/2 

n k = {g 2 k [ala k + (l + a k z )/2]+Al/A] . 

We will consider that the atoms are initially in an "X" -state, characterized by a (reduced) atom density operator 
whose non-zero elements are found only in the main diagonal and antidiagonal in the basis ||T) > |TI) , lit) > III)- This 
class of atom states has the property that the correspondingtwo-qubit density matrix preserves the "X" form when 
evolving under the action of certain system dynamics @,0,H, 0j- In this case [1,0], the reason is simple. When an 
atom transfers its excitation to the initially empty field, it leaves a signature in terms of the excitation in the field. 
Therefore, there cannot exist any coherence between the states ||T) and the states |||) , |||) unless such coherence 
existed initially. This is not the case, by definition, for the "X" -states, and therefore they will retain their "X" form 
under the assumed evolution whose form was motivated in the Introduction. 
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As pointed out in 5J, the "X" -class of states include the Bell states and the Werner states. Following 
focus on the Bell-like pure states, 



we will 



10(0)) = 



cos a |TT> - 
= cos a |T|) 



i/3 

4,3 'IT) 



s ") 'i' " III) , 



sm ae 



where < a < tt/2, < f3 < ir. (In Sec. |V]we will consider more general states.) As motivated above, we will assume 
that the initial state for the four qubit model is 



|*(0)) = M0))®|00) : 



|*(o)) = |v(o))®|oo) 



(4) 



respectively, where the abbreviated notation |0) ® |0) 6 = |00) has been used. Notice that Bell states can be recovered 
by setting a — ir/4 and (3 — 0, tt/2. The initial states (01 under the action of the operator JJ A ®Ub evolve as 



|*(t)) = x 1 |TTOO) + .t 2 |t|oi) + x 3 ||tio) 

+x 4 \H ll)+ar 5 |U00), 
|*(*)) = Vi IU 00) + y 2 |TT 00) 

+2/3 IU 10) + y 4 |U01) ! 



where the coefficients for the state ([5]) are given by 



%i = fA(t)fB(t) cos a, 

xi = fA(t)9B(t) cos a, 

X3 = 9A(t)fB(t) cos a, 

x 4 = g A (t)g B (t) cos a, 

x 5 = hA(t)hB(t)e %l3 sin a. 



The functions fk(t), gk(t), and hk(t) are given by 

/*(*) 



cos(rifet) - i— — sm(fifct) 



9k(t) = -i^e^'sin^t), 



h k (t) 



JA k t/2 



where the Rabi frequencies are flk = (gl + A^/4) 1 / 2 . 
Similarly, the state © will have the the coefficients 



Vi = fA(t)h B (t) cos a, 
2/2 = h A (t)f B (t)e 



i/3 



sin a, 



2/3 = gA(t)h B (t) cos a, 
yi = h A {t)gB(t)e tf3 sin a. 



(5) 
(6) 



(7) 
(8) 
(9) 
(10) 

(11) 



(12) 

(13) 
(14) 



(15) 
(16) 
(17) 
(18) 



III. ENTANGLEMENT DYNAMICS 



In this section we will analyze the entanglement evolution in the double JC model. As an entanglement measure 
we will use a wedge-product based measure introduced in [17| . which, for the two-qubit case, coincides with the well- 
known concurrence [15j |, and in the multiqubit case with the entanglement monotones [l8l |. This measure is defined 
for any number of subsystems, each having an arbitrary, but finite, dimension. 

By partitioning the total system into two we can compute the entanglement between these partitions. Consider 
some partition composed by P\ with dimension M and P 2 with dimension N (note that each partition could contain 
more than one physical subsystem). Assume a pure system defined by 

M N 

rn—l n—1 



4 



Cab< ^Aa-Bb C AB ,4£ Aa _ Bb 

1.0 K 7 r 1.0 k r, 

:\ / \ \ / \ 

0.8 A \ //'\\ 0.8 -W //V 

M // \\ - \\ // M 

0.6 - \\ /• \^ 0.6 - /; A 

M " A 0.4: >i " .< 

0.2- \\ // , X N °- 2 : \\ // \\ 

- ^ f v • ' " \/AxN// V^Va 



5 (a) 1 2 3 4 5 (b) 



FIG. 1: (a) The concurrence Cab between the atoms in state |\P(0)) as a function of the coupling parameter Q.At when 
gA/gB = 1 for a = 7r/4 (dashed line), a — n/6 (dash-dotted line); and AEAa-Bb for a = 7r/4 (solid line) and a = tv/6 (dotted 
line), (b) The concurrence when <7s/(?a = 2 for the parameter values a — ir/A (dashed line) and a — it/6 (dash-dotted line); 
and AEAa-Bb for a = 7r/4 (solid line) and a = 7r/6 (dotted line). 



where and {|n}} are orthonormal bases. In order to estimate the entanglement between partitions P\ and P2, 

we project onto the basis states of one of the partitions. To this end we define the unnormalized state 

\ipm) = (mltp). 

If the system can be written as a tensor product between a pure state in each partition, then all states \ip m ) are 
parallel. That is, \ip m ) = c m for all m = 1, . . .M, where c m denotes a c-number. If, on the other hand, the pure 
state (I19p is entangled, then at least two of the vectors, say \ip m ) and \if>i) are not parallel, and the degree to which 
they are not parallel is characterized by the "area" the vectors span. This area is given by the wedge product between 
the vectors, but as the wedge product, in general, is signed and complex, we take the absolute square of the area as 
a measure of the nonseparability between these two vectors. The square of the measure introduced in [17| can hence 
be written as the determinant 



A 2 {m,l) 



(tpm\tpm) (lpm\lpl) 
(4>l\4>m} (4>l\?Pl} 



Summing all contributions and using symmetry and the fact that the wedge product between a vector and itself 
vanish, the entanglement between P\ and Pi can finally be defined 



M M 



(20) 



m=l 1=1 



In the case of the double-JC model, the possible partitions are: (a) one qubit - three qubit partitions, A — Bab, 
B — Aab, a — ABb, and b — ABa; (b) two qubit - two qubit partitions, Aa — Bb, Ab — Ba, and AB — ab. 

In Figs. [T]and [2] we plot the evolution of the concurrence Cab between the two atoms A and B, and the evolution 
of AEAa-Bb (since < Cab < 1, we use 4EA a -Bb to scale it to essentially the same range of values as C may obtain) 
in the case of exact resonance (A^ = = 0). Note that while Cab represent only the remaining entanglement 
between the two atoms after the field states have been traced out, EAa-Bb represent the entire bipartite entanglement 
between the atom-field system Aa and Bb. 

First we consider different values of a. In Fig. [T] (a), corresponding to gA = Qb, we can observe that concurrence for 
the initial state |^(0)) evolve in a typical oscillating way between and 1 [I,!!,!!!], meanwhile EAa-Bb = sin 2 acos 2 a 
is invariant and equals 1/4. In Fig. [l](b) we show the time evolution of the concurrence between the atoms, and that 
of ^EAa-Bb-, for c/b/qa = 2. When gA ^ gs the concurrence is not evolving in a typical oscillatory manner as was 
pointed out in [12j. Meanwhile, EAa-Bb depends only on the initial state, namely on a, and not on the ratio gA/gB- 

Fig. [2]shows the evolution of the concurrence and AEAa-Bb for the state |3>(0)) when both cavity-atom systems are 
in exact resonance. In Fig. [2] (a) we can observe the so-called entanglement sudden death for different values of a. 
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FIG. 2: (a) The concurrence Cab between the atoms in state |$(0)) as a function of the coupling parameter Q^t when 
gA/gB = 1 for a = 7r/4 (dashed line), a = tv/6 (dash-dotted line); and 4EAa-Bb for a = 7r/4 (solid line) and a = tv/6 (dotted 
line), (b) The concurrence when <?s/(?a = 2 for the parameter values a — 7r/4 (dashed line) and a — tt/6 (dash-dotted line); 
and AEAa-Bb for a = 7r/4 (solid line) and a = tv/6 (dotted line). 



(In the next section we shall discuss this phenomenon in more detail.) In Fig|5] (b) we show the effect of a different 
ratio between 17,4 and gs- 

The difference between the evolution of the two states arises because state | \& (0) ) is evolving simultaneously in 
two closed manifolds, one consisting of the one-excitation, subsystem Aa manifold (consisting of states ||J, 00) and 
I J, J, 10)), the other consisting of the one-excitation, subsystem Bb manifold (consisting of states |J,| 00} an d ||| 01)). 
On the other hand, the state |3>(0)) is evolving only in one manifold, consisting of the states |T| 00), |f| 01)), 10), 
and (J,! 11). State |4>(0)) also has a ground-state component ||| 00) but this state does not evolve. 



IV. DISSIPATIVE DYNAMICS 



While the model presented above is closed and does not involve any dissipation, the atoms' evolution under dissi- 
pation can still be described. As can be seen from Eqs. (|7|)- (fl"0|) , (fl"2f . (fl~3|) . and (fl~5| - (fl~8|) the excitation of the atoms 
is transferred to the fields in a monotonic fashion during the time interval where < U,^t < it/2. In the resonant 
case (Aa — As = => ilk = gk), all the excitation will be transferred from the atoms to the cavity fields. Formally, 
during this time interval, one can then see the fields as a dissipativc channel for the atoms' excitation. On can 
subsequently map a dissipative evolution, e.g., spontaneous emission of the atom's excitation obeying an exponential 
decay oc exp(— Jkt'), onto the JC dynamics by the identification between the times t and t'\ exp(— jkt') — cos 2 (ilkt). 
Quite obviously, fifei — > ir/2 correspond to t' —> 00. Hence, if the entanglement between the atoms (after the fields 
are traced out) become zero in a time r < 7r/(2f2fe), then in the dissipative picture it vanishes in a finite time 
t' = 7 _1 ln(cos _2 [S!fcT]). This is indeed what happens for the state \<fi(0)) as it evolves. 

One may then ask why one state's entanglement vanishes in finite time while the others' does not. The reason is the 
fundamentally different way the states decay. The state | -0(0)) decays directly into the ground state |||). Whatever 
excitation is left in the atoms will still be in a superposition state, and such a statistical mixture between a Bell state 
and the ground state cannot be written as a convex sum of any separable states no matter to what extent the state 
has decayed. The state |</>(0)} decays to the ground state via the intermediate states |||) and ||f ). As the decay leaves 
different "signatures" in the reservoirs (the states |01) and |10), respectively), no coherence between these states is 
established. When the excitation of these intermediate states is large compared to the remaining coherence between 
the states |TT) and |||), the state can be written as a convex combination of separable states so the state is no longer 
entangled. This happens when 

tana < sin^S!^) = 1 - exp(-7fct') (21) 

as pointed out in 0, If cosa < sina, the atomic excitation is insufficient to excite the intermediate states |||) 
and ||t) to the extent that the entanglement between the atoms vanishes in a finite time. 

If the coupling constants gA and gs are different, say that gA > gs as in Fig. Q] (b) and [2] (b), the just presented 
dissipative picture is only valid as long as < Qa^ < tt/2, where Qk = gk when the atoms and cavities are resonant. 
For times longer than ir/(2gA) the excitation of atom A starts to revive again in the JC model, a phenomenon not 
compatible with dissipation. However, as seen from the figures, the behavior of the states for times fJ^i < ir/2 <^> 
£lst < 7r/4 « 0.79 is qualitatively the same as in the symmetric (gA = <?s) case. 



6 



V. AN ENTANGLEMENT INVARIANT 

The form of the chosen interaction Hamiltonian, which does not include any interaction between subsystems Aa 
and Bb, ensures that no entanglement is formed between these subsystems that was not already present in the initial 
state. Therefore, it is reasonable to expect an entanglement invariant to exist that measures the net entanglement 
between these subsystems. As FigQ]and Fig J5] suggest, we can find the invariant 

£ = E Aa _ Bb (22) 

valid for all parameter values, that is, even for nonresonant coupling and different atom-cavity coupling ratios. 

The introduced measure ([22]) does not depend on time for the initial pure "X" -states ([5]) and ©. The value of £ 
in both cases is 

£ = sin 2 a cos 2 a. (23) 
This value is proportional to the square of the two atoms' initial concurrence C flEl ]: 

4£ = C 2 . 

The result is expected, because as the Hamiltonian is chosen not to change the entanglement between Aa — Bb, it 
must remain equal to its initial value at all times. What is more significant is that for the initial generic state 

|&) = a |TT 00} + c 2 |1T 10) + c 3 |TI 01} + c 4 |U 11) + c 5 ||| 00} 

+di |TI 00} + d 2 ||t 00} + d 3 HI 10} + d 4 HI 01} , (24) 

we find that (|22p is still invariant during evolution under the action of V a ® Ub ■ Explicitly, 

£(!£<?)) = l c i c 4 - c 2 c 3 | 2 + \ad 3 - c 2 d 1 \ 2 + \c 3 d 3 - c 4 d 1 \ 2 

+ |cid 4 — c 3 d 2 \ 2 + \c2d4 — c 4 d 2 | 2 + \c1C5 — d\d 2 \ 2 

+ |c 2 c 5 - d 2 d 3 \ 2 + |c 3 c 5 - did 4 | 2 + |c 4 c 5 - d 3 d 4 \ 2 . (25) 

It is worth noticing that each cavity-atom in the state |£ 9 ) will evolve in the subspace {|f 0) , || 1) , || 0)} under the 
double JC Hamiltonian. Hence, the state can be seen as a two-qutrit, pure state. In this case equation (|2"5]) is 
just one ninth of the square of the two-qutrit concurrence introduced in [19J. 

It should also be noted that in the dissipative picture, the state |£ s ) models coupling to excited reservoirs, and in 
general the states cannot, at any time, be written as a product state between the atoms and the fields. Hence, some 
entanglement between atom and fields is already there at the start of the evolution. The states |$(0)) and |^(0)) are 
special cases of |£ g } where the atoms couple to initially empty reservoirs, a relevant but special case. 

When we study other partitions in (|20|1 . we have seen that for both states, any single term is zero only at discrete 
times. This means that at all times, except a discrete set of times of zero measure, all parts of the system become 
entangled in some degree through excitation transfer. This phenomenon is generic for entangled systems and has, 
e.g., been used to entangle subsystems that have never interacted through so-called entanglement swapping [20| . 
In a dissipative system, this entanglement spread is of course detrimental and may lead to complete elimination of 
entanglement. 

VI. DISCUSSION 

In this work we have discussed the entanglement dynamics for two excited atoms coupled to cavity field-modes 
through a Jaynes-Cummings Hamiltonian. We have discussed both the closed system dynamics and the dynamics 
if the cavities are viewed as reservoirs. We have discussed why different initial atomic states face different fates 
(asymptotic vs. sudden decay) with respect of their entanglement under dissipation. We have also shown that as 
expected, there exist an entanglement invariant valid for a large class of states, much larger than the class of so-called 
"X" -states in the closed system. When treating the fields as reservoirs (i.e., when tracing over the fields), some of 
the entanglement transferred to the cavities is ignored. The state </>(0)} transfers its excitation over a larger set of 
distinguishable field states (|01), |10), and |11)) than the state |^(0)) (that excites only the field states |01) and 1 10) ) , 
and therefore is is not so surprising that the former state may loose all of its entanglement through even a finite 
dissipation. 

Since the proposed (bipartite) invariant is a measure of the entanglement between the two systems Aa and Bb, and 
the measure is invariant to local unitary transformations, it will be a constant for any pure state, not only for the 
generic state (|2~4")l . but also for states with higher excitation. 
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